Extending Hudson's theorem to mixed quantum states 
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According to Hudson's theorem, any pure quantum state with a positive Wigner function is 
necessarily a Gaussian state. Here, we make a step towards the extension of this theorem to mixed 
quantum states by finding upper and lower bounds on the degree of non-Gaussianity of states with 
positive Wigner functions. The bounds are expressed in the form of parametric functions relating the 
degree of non-Gaussianity of a state, its purity, and the purity of the Gaussian state characterized 
by the same covariance matrix. Although our bounds are not tight, they permit us to visualize the 
set of states with positive Wigner functions. 



The Wigner representation of quantum states [l|, 
which is realized by joint quasi-probability distributions 
of canonically conjugate variables in phase space, has 
a specific property which differentiates it from a true 
probability distribution: it can attain negative values. 
Among pure states, it was proven by Hudson Q (and 
later generalized to multi-mode quantum systems by Soto 
et Claverie @) that the only states which have non- 
negative Wigner functions are Gaussian states [|[. The 
question that naturally arises 0] is whether this theorem 
can be extended to mixed states, among which not only 
Gaussian states may possess a positive Wigner function. 
A logical extension of the theorem would be a complete 
characterization of the convex set of states with posi- 
tive Wigner function. Although this question can be ap- 
proached by using the notion of Wigner spectrum [f| , a 
simple and operational extension of Hudson's theorem 
has not yet been achieved due to the mathematical com- 
plications which emerge when dealing with states with 
positive Wigner functions 

Motivated by the increasing interest for non-Gaussian 
states in continuous-variable quantum information the- 
ory (see, e.g., p) and the need for a better understand- 
ing of the de-Gaussification procedures for mixed states 
(see, e.g., Q), we attempt here an exploration of the set 
of states with positive Wigner functions using Gaussian 
states as a reference. More precisely, we consider the sub- 
set of such states that have the same covariance matrix 
as a reference Gaussian state. We obtain a partial solu- 
tion to the problem, by analytically deriving necessary 
conditions (bounds) on a measure of non-Gaussianity for 
a state to have a positive Wigner function. This set of 
conditions bounds a region in a three-dimensional space 
with coordinates being the purity of the state, the pu- 
rity of the corresponding Gaussian state, and the non- 
Gaussianity. As intuitively expected, the maximum de- 
gree of non-Gaussianity increases with a decrease in the 
purity of both the state and its Gaussian corresponding 
state. 

Before deriving the main results of this paper, let us 
recall a convenient representation of the trace of the prod- 
uct of two one- mode quantum states, p and p ', in terms 



of the Wigner representation |8j, 

Tv(pp') = 2nJ J dxdpW p (x,p) W p , (x,p) , (1) 

where W p is the Wigner function of the state p. For 
example, the purity of a state, p[p] = Tr (p 2 ), may be 
calculated with the help of this formula. For a state with 
a Gaussian Wigner function determined by the covari- 
ance matrix 7 and displacement vector d, the purity is 
simply p [pa] = (det 7) . The matrix elements of the 
covariance matrix of state p are defined as 



Hi = Tr({(r 4 - d,), (fj - dj)}p) 



(2) 



where f is the vector of quadrature operators f = (x,p) T , 
d = Tr(rp), and {■, ■} is the anticommutator. Note that 
we can put the displacement vector to zero with no loss 
of generality since the purity (and all quantities we will 
be interested in) does not depend on d. We will thus 
consider states centered on the origin in this paper. 

Our aim is to derive bounds on the non-Gaussianity, 
i.e., on the "distance" between a state p of purity p[p] 
possessing a positive Wigner function and the Gaus- 
sian state pa determined by the same covariance matrix. 
While there are different measures in the literature for 
quantifying the distance between two mixed states, we 
have chosen to use a recently proposed one @, 



$ [p, Pg] 



p[p] +p[pg] - 2Tr(pp G ) 
2p[p] 



(3) 



Although the quantity S[p, pa] is obviously not symmet- 
ric under the permutation of the two states, it is conve- 
nient for quantifying the non-Gaussian character of p in 
the sense that 8 £ [0, e], with e < 1, and 8 — is at- 
tained if and only ifp = pg- For one- mode states, it is 
conjectured in Ref. Q that e = 1/2. 

In a first step, we are going to derive bounds on the 
trace overlap Tr (ppg) for fixed values of p[pc] and p[p\. 
It will then be straightforward to express bounds on the 
non-Gaussianity 8 [p, pc] in terms of p[pc) and p[p] by 
using Eq. ©. 
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We use Eq.|T|) in order to reformulate the problem as 
an optimization problem that can be tackled with the 
method of Lagrange multipliers. More specifically, we 
need to extremize the functional /[Wp] = Tr (ppo) rep- 
resented by Eq. {I} with the constraint that the Gaus- 
sian Wigner function W PG and the positive function W p 
possess the same second moments. In order to simplify 
our derivation, we apply a symplectic transformation S 
on the states p and pa, giving SpS^ and SpcS^ respec- 
tively, in such a way that the Gaussian state becomes 
invariant under rotation in the x-p plane (i.e., becomes 
a thermal state). In this way the problem is reduced to 
a simpler but equivalent one, since the functional /[Wp] 
and the purities of the states remain invariant under S 
and since the positivity of W p is preserved. This last 
statement can be justified by the fact that the time evo- 
lution of a Wigner function under a quadratic Hamilto- 
nian can always be viewed as an affine transformation on 
the variables x andp Q. Furthermore, we claim that the 
function W p x which extremizes the functional I[W P ] is 
invariant as well under rotation in the x-p plane, and we 
will justify this assumption at the end of the derivation. 

After the application of the symplectic transformation 
and under the assumption of rotation-invariant solutions, 
the functions W p (r) and W PG (r) = 2^c e ~ r ^ 2C only de- 
pend on the squared radius r = x 2 + p 2 , and the func- 
tional /[Wp] is written in a simpler form as 

/>oo 

I[W P ] = Tr (pp G ) = 2ir 2 / W p (r) W PG (r) dr. (4) 
Jo 

The constrains that we impose on the function W p can 
be summarized as follows: 

(1) It is positive for the values of r belonging to some 
set s and zero elsewhere. 

(2) It is normalized, ir J s W p (r)dr = 1. 

(3) It has the same variance as the corresponding Gaus- 
sian state, pa 



(5) It is square integrable and continuous. 

This last requirement follows directly from the general 
property of Wigner functions, 



W(x,p)dp = (x\p\x) 



(7) 



Recall that a state can always be diagonalized in a basis 
of pure states, namely p = X^j^i \tpi) Since wave 

functions must satisfy the conditions of continuity and 



integrability in both position and momentum represen- 
tation, one concludes that a Wigner function of variables 
x and p, and more generally of any variable that is a con- 
tinuous function on these, e.g., r = x 2 +p 2 , has to satisfy 
the same requirements. 

Finally, let us stress that without the requirement of 
positive definitcness of the operator p, the set of condi- 
tions listed above is not sufficient to constrain the solu- 
tions Wp(r) to eligible Wigner functions. To our knowl- 
edge, there exists no operational criterion on phase-space 
functions ensuring that the operator p is physical (see 
[l(j ] for an extensive discussion). On the other hand, one 
can verify whether a quasi-probability distribution is un- 
physical by using a theorem which states that a square 
integrable and normalized function is an eligible Wigner 
function if its overlap with the Wigner function of every 
pure state is positive (Tl| . 

After having applied the method of Lagrange multipli- 
ers, we obtain the extremal solution 



W e p x (r) =A 1+ A 2 



1 



2ttC 



-r/2C 



+ A 3 r, 



(8) 



ir [w p (r)rdr=2C=l/(i[p G ]. (5) 

(4) It is such that the state p has purity p,[p], 

2tt 2 f W 2 (r)dr = n\p\. (6) p[p] ex - v\pg 



with the A J s being determined by conditions 2-4. Square 
integrability, condition 5, limits the class of possible func- 
tions Wp X {r) in Eq. ([5]) to those that have zero, one, or 
two positive roots denoted as tb (in the one- and two- 
root cases) and (in the two-root case). Furthermore 
the condition of continuity dictates that s = [r^rs] in 
the two-root case, s = [0, rs\ in the one-root case, and 
s = [0, oo] in the zero-root case. The latter case is the 
trivial one, where Wp X {r) coincides with W PG (r) and thus 
S [p, pa] vanishes. We treat the other two cases sepa- 
rately, and obtain two continuously connected branches 
of solutions for Wp x . The expressions that we obtain for 
Tr (ppa) ex and p\p] ex are highly non-linear, so that it is 
not possible to derive an analytic expression that directly 
connects the two quantities. Nevertheless, we are able to 
express the extremal solutions in the form of parametric 
functions. 

(I) Two roots : W^ x {r A ) = W^ x (r B ) = 0. We express 
the extremum purity /i[p] ea: and overlap Tr (ppc) ex in 
terms of the purity of the corresponding Gaussian state 
p\pc] and parameter a = (t\b — rA)p[pa], 



2 (a 2 -9 sinh(a)a + 2 (a 2 + 6) cosh(a) - 12) 



3a (a cosh (f) -2sinh(f))' 



^{pPgY x = PiPG^exp 

x (e a - 1) /a, 



a(a+e°(2a-3)+3) 
3(e°(a-2)+a+2) 



(9) 



(10) 



where < /k[/3g] £ 1- By imposing the condition rA > 
0, we obtain the bound < a < x r , with x r being the 
root of equation, 



e x (x-3) + 2x + 3 = 0. 



(11) 
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(II) One root : Wp X (rB) = 0. The extremal solution is defined by the following pair of parametric functions, 



4 ( e 2 ^(/? - 3) 2 + 8e^(/3 - 3) + /3(/3(2/3 + 9) + 12) - 9) 

u p = « pg o , 12) 

(2e/3(/3-3) + /3(/3 + 4) + 6) 2 

x 1 4(/3(cosh(/3) + 2)-3sinh(/3)) 

Tr( PPG ) = p[p G ] 2e/3(/3 _ 3)+/3(/3 + 4) + 6 , (13) 



where < /u[/?g] < 1 an d (3 — tbp[pg\- The range of 
the latter parameter is (3 > x r . 

We now need to show that, although we have only 
considered solutions Wp X with no angular dependence, 




FIG. 1: (Color Online) (a) Upper bound on non-Gaussianity 
5 [p, Pg] bx derived with the Lagrange multipliers method. The 
double thin orange line marks the boundary between the two 
branches of solutions. The dotted blue line indicates the in- 
tersection with the plane p, [p] = 1 (also noted 8 u ' ult in Fig. 
2) while the double thick blue line shows the intersection 
with the plane p. [pa] ~ 1- The red thick straight line de- 
notes the 'left' extremity of the surface, (b) Lower bound on 
non-Gaussianity 8 [p, pc] cs implied by the Cauchy-Schwarz 
inequality. We plot this lower bound only up to the inter- 
secting line of 5\p, pg\cs an( l 3[p,pG] ex (p[pa] = p [p] and 
5[p,p G ]=0). 



our result is general. If we waive this assumption and 
consider the most general case, we arrive to 

W?{x,p) = A 1+ A 2 ^e^ 2+ ^ c 

+A 3 x 2 + A 4 p 2 + A 5 xp, (14) 

which is the analog of Eq. ([5]) but allowing for an an- 
gular dependence. We can then apply a phase rotation 
on states p and pc in order to eliminate the term A 5 xp 
in Eq. (fT4|) . Such a rotation does not affect the cor- 
responding Gaussian state (since it is thermal) nor the 
trace overlap, so that the resulting extremal function be- 
comes symmetric by reflection with respect to the x or p 
axis in phase space. The conditions (x 2 ) = (p 2 ) = C for 
the function in Eq. (fT4]) is satisfied if A3 = A4,. Thus, the 
most general solution reduces to the rotation-invariant 
one, namely, Eq. ©. 

By using the derived bounds on the trace overlap 
[Eqs. Ij9"])-([T3"])], we plot in Fig. [ITa) the corresponding 
(upper) bounds on the non-Gaussianity S [p, pc\ x ■ By 
direct inspection, we conclude that the intersection of 
the plotted surface with the plane of pure states p,[p] = 1 
provides us with an upper bound on the non-Gaussianity 
of any state with a positive Wigner function and fixed 
covariance matrix (or, equivalently, fixed p[pc\) which 
is independent of its purity fi[p]. We denote it as the 
ultimate upper bound S u (p [pc]) and its parametric 
expression can be directly derived by setting p ex [p] = 1 
in Eqs. (O (USD ■ In Fig. [U we plot 5 uMt together with 
a lower estimation on this, S u , obtained by the con- 
vex combination of two symmetrically displaced coherent 
states. The tight ultimate upper bound on S [p, pa] must 
be located between these two curves. 

The bounds derived by using the Lagrange multipliers 
method confine S only from above. In order to obtain a 
lower bound on <5, we need to find an upper bound on the 
trace overlap. We achieve this by applying the Cauchy- 
Schwarz inequality on the Wigner representation of the 
trace overlap, [Eq. JIJ]. By using the definition of the 
purity, we arrive at 



Tr (pp G ) < v 7 P [pg] P [p] = Tr (ppg) c .s. ( 15 ) 

where CS stands for "Cauchy-Schwarz,". This bound, 
displayed in Fig. [IJb) , delimits together with the upper 
bound of Fig. [TJa), the region accessible for states with 
positive Wigner function in this 3D representation. Let 
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FIG. 2: (Color Online) Ultimate upper bound 5 u ' uU (blue 
solid line) on the non-Gaussianity of states with positive 
Wigner fuction as a function of the putity of corresponding 
Gaussian state fi[pa]. A lower estimate S Lult (black dashed 
line) on the ultimate upper bound that was obtained for a 
mixture of coherent states. The two curves limit the region 
where the tight ultimate upper bound on non-Gaussianity 
must be located. 

us note that this lower bound holds for states with both 
positive and negative parts Wigner functions. 

Let us now address the question of the physicality of 
the upper bound, namely the extremal solution Wp X . By 
resorting to Hudson's theorem, we can conclude that the 
intersections of the surface with the planes ji\p] = 1 and 
iApg] = 1 [ see single and double blue lines in Fig.QJa))] 
cannot correspond to physical states where 6^0. The 
only physical solution belonging to these lines is thus one 
point, namely, fi[p] — 1, p\pc] — 1, and 5 = 0. In order 
to test the physicality of the rest of the surface, we ap- 
plied the theorem mentioned above employing the eigen- 
states of the quantum harmonic oscillator as test pure 
states. From our analytical results on the first 40 num- 
ber states, we infer that the only functions Wf; x giving a 
positive overlap with every number state as n — > oo are 
the states with p[pc] — 0, that is, infinitely mixed states. 
Therefore, we conclude that the extremal solution of the 



form of Eq. (jSJ is unfortunately unphysical; hence, our 
bound is not tight. 

Finally, one may notice in Fig. H^a) that the left ex- 
tremity of the bound (red thick straight line) is on the 
left of the plane p[pc] — p[p\- The equation for this line 
can be easily derived, 

g 

p[p] = q^Pg} (16) 

and thus sets a lower bound on the purity of a mixed state 
given the purity of the corresponding Gaussian state. 
This bound has been derived in another context by Bas- 
tiaans [121 ] and has been proven to be the asymptotic 
form of an exact expression derived later by Dodonov 
and Man'ko (l3j in the context of purity bounded uncer- 
tainty relation. The exact bound is more strict than the 
bound in Eq . (116|) . and it is realized by positive Wigner 
functions [141 ] . This fact confirms again that our bound is 
unphysical but it also gives some evidence about the un- 
derlying link between Hudson's theorem and the Hcisen- 
berg uncertainty principle. 

In conclusion, we have found both upper and lower 
bounds on the non-Gaussianity of mixed states with pos- 
itive Wigner function. These bounds only depend on the 
purity and covariance matrix of these states, and an ul- 
timate upper bound can be derived that does not even 
depend on the purity, making it experimentally accessi- 
ble. An open question remains to derive tighter bounds 
for the non-Gaussianity. All our results apply to one 
single mode, so another natural question would be to in- 
vestigate the case of several modes. 
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